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PREFACE 
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The  manuscript  was  submitted  for  publication  on  September  17,  1980. 
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1.0  INTRODUCTION 

Studies  of  small  particles  entrained  in  flows  are  of  importance  in  a  number  of 
applications  at  the  Arnold  Engineering  Development  Center  (AEDC).  Both  velocity  and  size 
distributions  of  particulate  flows  are  of  interest.  The  goal  of  the  research  reported  here  was 
the  development  of  a  particle-sizing  interferometer  (PSI)  capable  of  providing  both  velocity 
and  size  information  for  use  in  AEDC  facilities. 

Specifically,  a  requirement  existed  for  particle  size  instrumentation  for  the  Engine  Test 
Facility  (ETF)  subscale  icing  facility,  RID.  Conditions  of  flow  velocity  (about  100  meters 
per  second),  particle  number  density  (a  few  hundred  per  cubic  centimeter),  and  range  of 
particle  diameters  (5  to  80  micrometers)  were  appropriate  for  interferometric  sizing.  A  PSI 
was  developed  for  these  conditions,  and  this  report  documents  the  development. 

Two  different  systems  were  evaluated,  one  using  on-axis  light  collection  and  an  analog 
signal  processor,  and  the  other  using  off-axis  light  collection  and  electronics  based  on  a 
commercially  available  high-speed  transient  signal  digitizer.  The  off-axis  light  collection 
arrangement  was  designed  to  provide  a  smaller  observation  volume  when  it  was  suspected, 
during  evaluation  of  the  on-axis  system,  that  signals  were  being  received  from  more  than  one 
particle  at  a  time. 

The  optics  and  electronics  of  the  two  systems  are  described  in  Sections  3.0  and  4.0  of  this 
report,  following  a  brief  review  of  the  theory  of  interferometric  sizing  in  Section  2.0. 
Corrections  to  the  raw  data  to  eliminate  system  bias  toward  certain  size  ranges  are  the 
subject  of  Section  5.0.  The  system  monosize  response  and  its  implications  are  discussed  in 
Section  6.0.  Laboratory  evaluation  and  wind  tunnel  data  results  are  presented  in  Section  7.0, 
and  the  performance  capabilities  of  interferometric  particle  sizing  are  discussed  with  the 

conclusions  of  Section  8.0. 

2.0  THEORY  OF  INTERFEROMETRIC  PARTICLE  SIZING 

The  theory  of  interferometric  particle  sizing  has  been  developed  and  reported  by  a 
number  of  workers  (Refs.  1  through  8).  The  optical  system  used  is  similar  to  a  high  quality, 
cross-beam  laser  velocimeter  and  is  outlined  in  Fig.  1.  A  laser  beam  is  split  into  two  beams 
either  by  glass  blocks  or  by  a  Bragg  cell.  The  two  beams  are  crossed  and  focused  at  the  same 
location  to  create  planar,  parallel  interference  fringes.  When  a  single  particle  crosses  the 
interference  fringes,  the  scattered  light  is  collected  by  a  lens  and  transmitted  to  a 
photodetector,  where  it  is  converted  to  an  electronic  signal  for  processing.  A  pinhole  or  slit 
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placed  in  front  of  the  photodetector,  through  which  the  collected  light  is  focused,  eliminates 
stray  signals  originating  outside  the  beam  crossover  region.  The  collection  lens  may  be 
located  either  on  the  optical  axis,  as  shown  in  Fig.  1,  or  at  an  angle  off-axis.  For  on-axis  light 
collection,  some  arrangement  of  beam  stops  or  mask  must  be  employed  on  the  collection 
lens  to  prevent  the  laser  source  beams  from  entering  the  photodetector.  A  typical  scattered 
light  waveform  is  also  shown  in  Fig.  1.  The  waveform  may  be  thought  of  as  having  two 
components:  a  mean,  or  “pedestal”  component  (shown  as  a  dashed  line),  to  which  is  added 
an  oscillating,  or  “a-c”  component.  The  “visibility”  of  the  fringes  in  the  waveform  is 
defined  as  the  ratio  of  the  magnitude  of  the  a-c  component  at  any  point  along  the  waveform 
to  the  pedestal  value  at  that  same  point.  The  visibility  may  vary  from  point  to  point  along 
the  waveform. 


a-c  Magnitude 
Pedestal 


Fringes 


Bragg 


Pinhole 


Electronics 


Figure  1.  Basic  schematic  of  optical  system. 

i 

The  scattered  light  signal  output  by  the  photodetector  may  be  approximated  in 
mathematical  form  as 


Jci  - 


ip0  (D)  exp  [-2/b2  (x2  +  7 2  +  z2  a2/4)] 
x  [cosh  (2  y z  a/b2)  T  Vo  cos (^2iry/ 8  2tt  fB 


(1) 


where  x,y,  and  z  are  coordinates  of  the  particle  location  with  respect  to  the  axes  shown  in 
Fig.  1,  b0  is  the  laser  beam  waist  radius  in  the  crossover  region,  a  is  the  full  beam  separation 
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angle  (assumed  small),  fa  is  the  Bragg  cell  driving  frequency  (set  equal  to  zero  if  a  glass  block 
beam  splitter  rather  than  a  Bragg  cell  is  used),  6  is  the  fringe  spacing,  VD  is  the  maximum,  or 
peak  value  reached  by  the  visibility  along  the  waveform,  and  t  is  time.  The  variable  ipo(D)  is 
the  value  of  the  pedestal  when  the  particle  is  centered  at  x  =  y  =  z  =  0.  It  is  a  function  of 
the  particle  diameter,  D,  laser  power,  light  collection  optics  geometry,  and  photodetector 
and  electronic  gain;  but  for  the  interest  here,  only  its  dependence  on  particle  diameter,  D,  is 
indicated . 

The  ratio  of  the  a-c  (i.e.,  the  cosine)  magnitude  to  the  pedestal  gives  the  visibility,  V,  as 
follows: 


V  =  Vo/cosh(2y  za/b2)  (2) 

A  previously  published  technical  report  (Ref.  9)  has  shown  how  the  peak  visibility,  VD, 
may  be  related  to  the  particle  size  for  a  given  collection  optics  geometry.  The  theory  was 
developed  using  Fraunhofer  diffraction  and  therefore  applies  to  light  collected  on-axis  in  the 
forward  direction.  The  size  of  the  collection  aperture  and  beam  stop  was  taken  into  account. 
Other  authors  (Refs.  4  through  7)  have  used  the  Mie  theory  of  scattering  from  spheres  to 
predict  visibility  versus  particle  diameter  for  off-axis  light  collection  and  for  particles  too 
small  for  scalar  diffraction  theory  to  apply.  A  recently  published  paper  (Ref.  10)  reports  the 
use  of  ray  optics  to  relate  off-axis  visibility  to  droplet  diameter  for  large  refracting  droplets. 

Curves  relating  peak  visibility  to  the  ratio  of  particle  diameter  to  fringe  spacing  (D/6)  are 
shown  in  Fig.  2  for  two  different,  on-axis  receiving  apertures.  The  dashed  curve  to  the  left  is 
for  the  theoretical  case  where  all  of  the  forward-scattered  light  is  collected  (corresponding  to 
an  infinite  collection  lens  diameter  and  no  beam  stops).  The  solid  curve  to  the  right  is  for  the 
aperture  created  by  placing  on  the  collection  lens  a  specialized  mask  which  admits  light  only 
through  two  60-deg,  wedge-shaped  openings  at  the  top  and  bottom  of  the  lens  as  sketched  in 
the  upper  right  of  the  figure.  This  curve  was  calculated  using  scalar  diffraction  theory  [Eq. 
(13),  Ref.  9]  and  assumes  spherical  scatterers  with  diameters  much  larger  than  the 
illuminating  wavelength.  Each  of  the  two  curves  shown  is  useful  for  sizing  over  roughly  a 
10:1  size  range.  At  the  upper  end  of  the  size  range,  ambiguity  results  from  the  fact  that  more 
than  one  particle  size  has  the  same  visibility.  At  the  lower  end  of  the  size  range,  the  curve 
becomes  too  flat  and  insensitive  to  size  changes  to  be  usable. 

The  fringe  spacing  is  chosen  for  the  particular  particle  size  range  of  interest  and  the 
applicable  visibility  curve.  For  example,  when  a  60-deg  mask  is  used,  the  value  D/6  =  2  lies 
at  the  upper  end  of  the  usable  size  range,  whereas  for  the  dashed  curve  of  Fig.  2,  the  value 
D/5  =  1  lies  near  the  upper  end  of  the  size  range.  Therefore,  to  cover  a  given  size  range,  the 
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60-deg  mask  uses  a  fringe  spacing  about  half  that  of  a  system  having  a  visibility  curve  similar 
to  the  dashed  one. 


D/6 

Figure  2.  Theoretical  and  experimental  comparison,  visibility 
versus  D/6  for  on-axis  light  collection. 

Experimental  data  obtained  using  glass  spheres  placed  on  a  glass  slide  are  shown  for 
comparison  with  the  60-deg  mask  curve.  Each  glass  sphere  was  sized  using  the  distance  from 
the  center  to  one  of  the  concentric  rings  in  its  single  beam  diffraction  pattern  for  comparison 
with  the  corresponding  visibility  measurement. 

For  light  scattered  near  the  optical  axis  in  the  forward  direction  from  large  particles  such 
as  those  in  the  size  range  of  interest  here,  there  is  minimal  dependence  on  the  scattering 
particle’s  index  of  refraction.  Experimental  data  agree  well  with  the  results  of  scalar 
diffraction  theory,  arrived  at  independently  of  the  refractive  index. 

A  curve  relating  visibility  to  droplet  diameter  for  a  case  where  the  scattered  light  is 
collected  off-axis  is  shown  in  Fig.  3.  The  collection  aperture  used  is  shown  in  the  upper  right 
of  the  figure.  The  collection  lens  is  centered  12  deg  off-axis  in  a  direction  perpendicular  to 
the  plane  of  the  beams;  the  lens  diameter  subtends  an  angle  of  1 1.4  deg  at  the  center  of  the 
beam  crossover  region.  The  curve  was  calculated  using  Mie  theory  for  the  case  of  water 
droplets  (refractive  index  =  1.33),  an  illuminating  wavelength  of  X  =  0.5145  fim,  and  a 
beam  separation  angle  of  or  =  2.7  deg,  which  gives  a  fringe  spacing  of  6  =  11  fim.  The 


8 


AEDC-TR-B0-50 


illuminating  beams  were  assumed  to  be  linearly  polarized  in  the  x  direction  (i.e.} 
perpendicular  to  the  plane  of  the  beams).  A  curve  calculated  for  the  case  of  oleic  acid 
droplets  (refractive  index  =  1.46)  closely  follows  the  curve  for  water  in  the  range  of 
diameters  from  5  to  30  fim. 


Figure  3.  Theoretical  and  experimental  comparison,  visibility 
versus  diameter  for  off-axis  light  collection. 


The  experimental  data  shown  were  obtained  using  water  droplets  and  oleic  acid  droplets 
produced  by  a  Berglund-Liu  vibrating  orifice  particle  generator.  The  diameter  of  the 
droplets  in  the  monosize  (monodisperse)  stream  produced  by  the  generator  could  be 
calculated  accurately  from  the  liquid  flow  rate  and  the  vibration  frequency  of  the  orifice 
(Ref.  11).  The  data  points  shown  are  centers  of  distribution  (histogram  modes)  measured  in 
response  to  monosize  droplets  flowing  through  the  beam  crossover  region. 

The  experimental  data  fall  noticeably  below  the  theoretical  curve  for  the  smaller  droplet 
diameters.  This  is  probably  due  to  imperfect  fringe  contrast  in  the  beam  crossover  region. 
The  calibration  curve  which  was  used  in  wind  tunnel  tests  to  relate  visibility  to  particle  size 
for  the  off-axis  system  is  shown  as  the  dashed  straight  line  and  is  based  on  the  experimental 
data. 
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3.0  OPTICS 


3.1  GENERAL 

Both  an  optical  system  using  on-axis  light  collection  and  an  optical  system  using  off-axis 
light  collection  were  built  and  evaluated.  Both  systems  consisted  of  two  major  units,  one  a 
transmitting  optics  package  for  generating  a  set  of  high-contrast  fringes,  and  the  other  a 
collection  optics  package  for  collecting  the  scattered  light.  In  both  systems  the  laser  beams 
were  crossed  and  focused  at  a  point  100  cm  in  front  of  the  final  transmitting  lens  to  form  the 
interference  fringe  set.  The  collection  lenses  were  a  set  of  15-cm-diam,  f5  lenses,  giving  a 
total  system  “throw”  of  1.75  m  between  the  transmitting  and  receiving  units. 

The  on-axis  system  used  a  15  mw,  helium-neon  laser  operating  in  the  TEM(K1  mode  at  a 
wavelength  of  0.6328  pim,  a  Bragg  cell  driven  at  15  MHz  to  split  the  beams,  and  a  system  of 
lenses  to  cross  and  focus  the  beams  at  a  common  point  to  produce  a  set  of  moving  fringes.  In 
the  wind  tunnel  application  the  fringes  moved  in  a  direction  opposite  to  the  direction  of 
particle  flow  to  accommodate  the  electronics.  In  the  laboratory,  a  stationary  particle  could 
be  positioned  at  the  crosspoint  with  the  fringes  moving  past  the  particle  to  produce  signals 
for  calibration  purposes.  The  primary  fringe  spacing  was  designed  for  26  ptm  with  the 
capability  of  obtaining  other  selected  fringe  spacings  by  changing  lenses.  The  beam  waist 
diameter  was  150  pun. 

Collection  optics  for  the  on-axis  system  included  the  specialized  mask  with  60-deg  wedge- 
shaped  openings,  a  200-ptm-diam  pinhole  to  block  light  originating  at  points  other  than  the 
beam  crossover  region,  and  an  RCA  8644  photomultiplier  tube. 

The  off-axis  system  used  a  150-mw  Argon  ion  laser  operating  at  a  wavelength  of  0.5145 
ptm.  The  beams  were  split  using  path-length-compensated  glass  blocks,  which  produced 
stationary  fringes  in  the  beam  crossover  region.  The  fringe  spacing  used  was  1 1  /im,  and  the 
beam  waist  diameter  was  200  pun. 

Scattered  light  for  the  off-axis  system  was  collected  in  accordance  with  the  arrangement 
of  Fig.  3.  The  light  was  focused  through  a  slit  120  ptm  wide  before  being  detected  by  an  RCA 
C30872  avalanche  photodiode. 

3.2  OBSERVATION  VOLUME 

Interferometric  sizing  requires  that  there  be  only  one  particle  at  a  time  in  the  field  of  view 
of  the  optical  system.  Extraneous  signals  from  large  particles  crossing  the  beams,  even 
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outside  the  beam  crossover  region,  can  interfere  with  the  measurement  of  a  small  particle 
located  within  the  crossover  region.  This  leads  to  the  concept  of  a  loosely  defined 
‘/observation  volume”  which  must  be  empty  of  additional  particles  when  a  particle  is  in  the 
crossover  region  and  a  measurement  is  being  made.  For  a  first-order  approximation  of  the 
requirements  which  must  be  imposed  on  the  size  of  the  observation  volume  for  a  given 
density,  n,  of  particles  in  the  flow,  the  assumption  is  made  that  the  observation  volume,  Vob, 
is  cylindrical,  of  radius  ba  equal  to  the  beam  radius  in  the  crossover  region,  and  of  length 
zob  ,  so  that  VDb  =  irb£  z<,b.  It  is  further  assumed  that  the  probability  of  occurrence  of 
multiple  particles  is  sufficiently  small  if  the  expected  number  of  particles  at  any  time  in  the 
observation  volume,  given  by  nVob,  is  0.1  or  less.  (Using  a  Poisson  distribution  of  events  to 
calculate  the  probability  of  finding  two  or  more  particles  in  the  observation  volume,  given 
that  there  is  already  at  least  one  there  to  provide  a  signal,  gives  S  percent  for  an  expected 
number  of  0. 1 ,  1  percent  for  an  expected  number  of  0.02,  and  42  percent  for  an  expected 
number  of  1.0.)  Then,  given  a  particle  density  in  the  flow  of,  say,  200  particles/cm3,  an 
upper  limit  may  be  set  for  the  size  of  the  volume,  as  follows: 

Vob  <  0.1/n  =  0.1/200  cm3  ^ 

<  0.5  mm® 

3.2.1  On-Axis 

Reducing  the  observation  volume  to  this  size  for  the  on-axis  system  required  design  of 
the  specialized  lens  mask  referred  to  above.  As  noted,  to  cover  a  given  range  of  particle 
diameters,  a  smaller  fringe  spacing  is  used  with  the  60-deg  mask  than  would  be  used  with  a 
system  having  a  curve  similar  to  the  dashed  curve  shown  in  Fig.  2.  This  is  an  advantage  in 
that  it  permits  one  to  use  a  smaller  focused  beam  radius  while  retaining  several  fringes  within 
the  region  of  beam  overlap.  The  observation  volume  is  correspondingly  reduced,  being 
proportional  to  the  beam  radius  squared. 

A  second  advantage  of  the  60-deg  mask  is  a  reduction  in  dynamic  range  of  signal 
magnitudes,  which  not  only  can  relax  requirements  on  the  electronics  but  also  makes  it  easier 
for  the  photodetector  pinhole  to  discriminate  against  extraneous  large  particle  signals,  thus 
shortening  the  observation  volume.  Plots  of  the  pedestal  ipo(D)  and  a-c  magnitude,  V0ip0(D), 
for  centered  particles  versus  D/6  are  shown  in  Fig.  4  for  this  mask.  These  curves  were 
calculated  during  the  course  of  evaluating  Eq.  (13)  of  Ref.  9  for  the  visibility.  For 
comparison,  a  dashed  line  shows  the  D2  dependence  of  the  pedestal  that  would  be  obtained 
if  all  of  the  diffracted  light  were  collected.  The  two  forward  lobes  (Airy  disks)  of  the 
diffraction  pattern  (one  from  each  illuminating  beam),  which  contain  most  of  the  diffracted 
light  power  and  cover  a  decreasing  angular  region  in  the  forward  direction  with  increasing 
particle  diameter,  are  blocked  by  the  mask  for  larger  particles,  giving  the  attenuated  curves 
of  Fig.  4. 
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Note  that  the  curves  of  Fig.  4  may  be  experimentally  determined  by  placing  particles  in 
the  center  of  the  beam  crossover  region  and  measuring  the  response  caused  by  the  moving 
fringes.  The  a-c  and  pedestal  components  are  separated  in  the  electronics  by  filtering  and 
may  be  measured  separately.  For  a  specific  system  (including  specific  laser  power,  photo¬ 
detector  gain,  and  electronic  processor  gain)  the  vertical  axis  may  be  given  specific  values. 
An  electronic  threshold,  or  “gate”  current  is  selected  as  shown  by  the  example  in  Fig.  4.  The 
ratio  of  signal-to-threshold  current  is  then  used  to  determine  the  “sample  volume”  cross- 
sectional  area  as  a  function  of  particle  size,  as  will  be  discussed  later. 


Figure  4.  Pedestal  and  a-c  magnitude  versus  D/5  for  on-axis 
light  collection. 

The  length,  zob,  of  the  observation  volume  is  determined  by  the  aperturing  effect  of  the 
60-deg  mask  in  combination  with  the  photodetector  pinhole,  an  effect  complicated  by  the 
varying  angular  distribution  of  the  scattered  light  for  different  particle  sizes.  A  limitation  on 
how  small  the  pinhole  may  be  made  in  order  to  reduce  Zob  is  that  it  is  desirable  that  the 
pinhole  not  block  any  light  scattered  from  the  particle  being  sized,  since  significant  blocking 
of  the  diffraction  pattern  beyond  the  lens  mask  alters  the  visibility  from  that  shown  by  the 
curve  of  Fig.  2.  A  pinhole  diameter  of  200  pm  was  selected,  which  with  this  system  gives  a 
length,  zob,  experimentally  determined  to  be  about  20  mm.  Since  bc  =  0.075  mm,  the 
observation  volume  of  the  system  is  about 
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V°b  =  *bo  7'ob 

-  f  (0.075) 2  (20)  =  0.35  mm3  (4) 

which  is  within  the  requirement  set  by  Eq.  (3). 

3.2.2  Off-Axis 

The  arrangement  of  collection  lens,  horizontal  slit,  and  photodiode  aligned  to  collect 
light  from  an  off-axis  angle  greatly  reduces  the  length,  Zob,  of  the  observation  volume.  This 
is  an  important  advantage  of  an  off-axis  arrangement.  Figure  5  shows  an  experimental  plot 
of  collected  light  scattered  from  a  glass  bead  versus  the  z  coordinate  of  the  bead  for  the  off- 
axis  arrangement.  The  120-/*m  slit.was  in  place  for  the  experiment.  A  spherical  glass  bead  of 
45-jim  diameter  was  placed  on  a  glass  slide  and  moved  by  a  micrometer  positioner  in  the  z 
direction.  The  coordinate  z  =  0  represents  the  center  of  the  beam  crossover  region,  although 

o  Experiment 
’ - Gaussian  Curve  for 


-1.5  -1.0  -0.5  •  0  0.5  1.0  1.5 


Coordinate  z.  mm 

Figure  5.  Collected  light  versus  particle  z  coordinate, 
off-axis  light  collection. 

the  bead  was  illuminated  by  only  one  beam  for  this  measurement.  The  collected  light  curve  is 
approximately  Gaussian;  for  comparison,  a  Gaussian  curve  with  standard  deviation  of  0.412 
mm  is  shown  as  a  dotted  line.  Very  similar  plots  were  obtained  for  beads  of  40-  and  30-yxm 
diameter.  From  the  plot  shown,  Zob  for  this  particular  off-axis  case  is  approximately  2.5  mm, 
giving  an  observation  volume 
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Voh  =  nhl  Zob 

=  it  (0.1)2  (2.5)  =  0.079mm3  (5) 

which  is  more  than  four  times  smaller  than  that  for  the  on-axis  case.  It  could  have  been  made 
yet  smaller  had  the  beam  diameter  been  reduced  to  that  of  the  on-axis  system. 

A  disadvantage  of  the  off-axis  arrangement  is  that  the  mechanism  by  which  the  light  is 
cut  off  for  particles  passing  at  z  coordinates  other  than  z  =  0  is  a  partial  or  total  blocking  of 
the  scattered  light  by  the  slit  in  front  of  the  photodetector.  Partial  blocking  of  the  light  alters 
the  visibility  from  that  shown  by  the  curve  in  Fig.  3.  The  result  seen  when  droplets  are 
allowed  to  freely  flow  through  the  crossover  region,  passing  at  random  z  coordinates,  is  a 
broadened  instrument  response  (i.e.,  a  recording  of  a  spread  of  particle  sizes  centered  on  the 
correct  value  when  a  monosize  stream  of  droplets  is  introduced). 

A  curve  showing  the  pedestal  ipo  for  centered  droplets  versus  the  droplet  diameter  is 
shown  in  Fig.  6  for  the  off-axis  arrangement.  This  curve  was  calculated  using  Mie  theory,  a 


Figure  6.  Pedestal  versus  diameter  for  off-axis 
light  collection. 
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droplet  refractive  index  of  1.33,  X  =  0.5145  jtm,  and  the  illumination  polarized  in  the  x 
direction. 


4.0  ELECTRONICS 

Different  electronic  signal  processing  and  data  acquisition  schemes  were  used  with  the 
two  optical  systems.  An  analog  signal  processor  was  used  with  the  moving  fringes  of  the  on- 
axis  optics,  and  a  fully  digital  approach  was  used  with  the  stationary  fringes  of  the  off-axis 
system. 

4.1  ANALOG  SYSTEM 

An  analog  electronic  signal  processor  for  measuring  the  visibility  of  the  scattered  light 
waveforms  was  built  on  the  basis  of  a  modification  of  the  block  diagram  of  Fig.  13  of  Ref.  9. 
The  modification  was  required  because  particle  velocities  would  have  created  frequencies 
beyond  the  bandwidths  of  existing  analog  dividers.  In  the  Final  design,  the  photodetector 
signal  was  amplified  and  sent  simultaneously  to  a  bandpass  filter  and  to  a  low-pass  filter  for 
separation  of  the  high  frequency  a-c  component  from  the  low  frequency  pedestal.  The  a-c 
component  was  envelope-detected  to  produce  a  signal  proportional  to  the  a-c  magnitude. 
The  a-c  and  pedestal  components  were  then  both  peak-detected  to  store  the  maximum  values 
reached  individually  by  the  a-c  and  pedestal  components  during  waveform  passage.  The 
ratio  of  these  two  values  was  measured  by  an  analog  divider  and  then  converted  to  digital 
form  by  an  analog-to-digital  converter.  This  ratio  is  a  valid  measure  of  the  peak  visibility, 
V0,  only  for  particles  passing  near  the  center  of  the  beam  crossover  region.  The  consequence 
is  that  a  high  electronic  threshold  must  be  used  to  ensure  that  only  particles  passing  near 
center  are  sized,  which  in  turn  means  that  smaller  particles  do  not  get  sized  because  they  do 
not  scatter  sufficient  light  for  their  signals  to  be  above  the  high  threshold.  The  result  is  a 
reduced  range  of  particle  sizes  that  can  be  covered  with  a  single  fringe  spacing.  This 
limitation  is  avoided  in  the  off-axis  system,  where  the  optics  are  designed  to  collect  light  only 
from  a  narrow  region  near  the  center  of  beam  crossover. 

A  microcomputer  based  on  the  8080  microprocessor  was  developed  to  take  the  individual 
visibility  readings  output  by  the  analog-to-digital  (A/D)  converter  and  put  them  in 
histogram  form.  Five  binary  bits  were  used  at  the  A/D  converter  output,  giving  thirty-two 
different  possible  visibility  readings,  numbered  0  to  31  and  corresponding  to  equal 
increments  of  visibility  from  0  to  1 .0.  Each  of  the  32  resolvable  visibility  readings  was  used 
as  the  address  of  a  memory  location  in  the  microcomputer.  The  content  of  each  memory 
location  was  set  to  zero  at  the  beginning  of  the  data  run.  Then  as  each  visibility  reading,  or 
address,  was  accepted  by  the  microprocessor,  the  content  of  the  corresponding  memory 
location  was  incremented  by  one  count.  At  the  end  of  the  data  run,  each  of  the  memory 
addresses  then  contained  the  raw  count  for  one  bin  of  the  histogram. 
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A  separate,  tabletop  computer  with  graphics  display  was  used  to  control  data 
acquisition,  including  control  of  a  counter  with  internal  clock  to  measure  the  total  time 
required  for  a  data  run.  The  content  of  each  of  the  32  memory  addresses  in  the 
microcomputer,  as  well  as  the  counter  reading,  was  accepted  by  the  computer  at  the  end  of 
the  data  run  via  IEEE-488  standard  interface  bus  cable  connections.  The  computer  then 
made  corrections  to  the  raw-count  histogram  as  will  be  described  in  Section  5.0  and 
displayed  the  corrected  size  distribution  histogram  on  the  graphics  screen.  Raw  data  could 
also  be  stored  on  magnetic  tape  for  later  analysis. 

4.2  DIGITAL  SYSTEM 

if 

The  signal  output, by  the  photodiode  in  the  off-axis  system  was  amplified  using  the  plug¬ 
in  vertical  amplifier  of  .an  oscilloscope.  The  waveform  was  then  stored  and  digitized  by  a 
commercially  available  high-speed  transient  digitizer  for  subsequent  software  analysis  using 
a  PDP-11  computer. 

The  computer  program  which  was  used  calculates  the  visibility  at  the  peak  of  the 
waveform  from  the  equation 


a-c  magnitude 
pedesLal 


(I 


->/* 


(I  +  I  )/2 

m  x  mn  / 


(6) 


where  1^  is  the  maximum  signal  current  occurring  at  the  peak  of  a  bright  fringe  and  I,™  is 
the  average  value  of  the  two  adjoining  dark  fringes.  The  waveform  is  rejected  for 
measurement  if  its  maximum  value  exceeds  the  limits  of  the  digitizer’s  storage  array,  or  if  the 
value  of  the  pedestal  (Imx  +  Imn)/2  is  not  greater  than  a  preselected  threshold  value.  The 
threshold  value  is  chosen  to  provide  adequate  resolution  of  measurement  for  small  signals. 
The  visibility  readings  are  sorted  into  bin  counts  by  the  program  for  display  in  histogram 
form  on  a  graphics  terminal  in  a  manner  similar  to  that  of  the  on-axis  electronics  display. 
The  histograms  could  also  be  stored  on  magnetic  disk  for  later  analysis.  A  commercial 
electronic  counter  was  used  to  measure  the  run  time. 


5.0  DATA  CORRECTIONS 


5.1  GENERAL 

Particle  size  distribution  data  are  often  presented  in  histogram  form.  A  hypothetical 
distribution  is  shown  by  the  curve  of  Fig.  7  with  some  of  the  histogram  bars  which  might  be 
used  to  represent  the  distribution.  The  distribution  curve  f(D)  is  a  plot  of  dn/dD,  where  n  is 
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Figure  7.  Example  log-normal  particle  size  distribution. 


the  total  number  density  of  particles  in  the  flow  in  particles/cm3  and  D  is  the  particle 
diameter  in  micrometers.  Units  of  f(D)  are  therefore  particles/cm3/^m.  The  area  under  the 
curve  is  equal  to  the  total  number  density  as  follows: 


f  CD)  dD 


n  particles/cm'3 


Sometimes  the  size  distribution  curve  is  normalized  so  that 


(7) 


f  (D)dD  =  1 


In  discrete,  or  histogram  form,  dn/dD  becomes  An/AD  and 


AD  =  £A" 


n  particlcs/cm^ 


(8) 


(9) 


The  output  of  either  data  acquisition  system  described  in  this  report  is  a  particle  size 
distribution  presented  in  histogram  form  on  a  graphics  terminal.  Certain  corrections  must  be 
made  to  the  uncorrected,  or  raw  histogram  bar  heights  before  the  histogram  is  displayed,  in 
order  that  it  be  in  the  proper  form,  An/AD.  Each  uncorrected  bar  height  represents  simply 
the  number  of  particles  counted  having  visibilities  falling  within  the  visibility  range 
represented  by  the  bar,  or  “bin”  width.  Each  raw  count  must  be  divided  by  the  volume  of 
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flow  sampled  to  obtain  An  particles/ cm3  and  then  divided  by  the  bin  width,  AD,  to  obtain 
An/AD. 

In  the  electronics  used,  each  histogram  bin  has  equal  width  on  the  visibility  axis. 
Reference  to  Fig.  2  shows  that  the  visibility  curve  for  the  on-axis  case  is  not  linear  and 
therefore  that  bin  widths  will  not  be  equal  on  the  size  axis. 

S.2  SAMPLE  VOLUME 

The  volume  of  the  flow  from  which  the  size  distribution  is  taken  is  called  the  sample 
volume.  It  is  a  cylinder  of  length  vtr  and  base  area  A,  where  v  is  the  flow  velocity,  tr  is  the 
time  during  which  the  electronic  processor  is  sampling  the  flow,  and  A  is  the  cross-sectional 
area  normal  to  the  flow  over  which  the  processor  actively  responds  to  measure  signals. 

The  active  sample  time,  tr,  is  given  by  the  equation 


where  T  is  the  total  elapsed  time  from  the  beginning  of  the  run  until  the  last  particle 
measurement  in  the  sample  is  made,  Nt  is  the  total  number  of  measurements  made,  and  t  is 
the  processor  downtime  between  measurements  (i.e.,  the  time  after  each  measurement 
during  which  the  processor  is  off-line  to  record  the  data  point).  If  a  particle  passes  through 
the  crossover  region  during  the  downtime,  the  processor  does  not  respond  to  it,  and 
therefore,  the  time  N,r  must  be  subtracted  from  T  in  the  determination  of  the  volume  of  air 
sampled. 

The  active  cross-sectional  area,  A,  varies  with  particle  size.  Although  small  particles 
must  pass  close  to  the  center  of  the  beam  crossover  region  for  the  scattered  light  signal  to 
exceed  the  processor’s  threshold  and  trigger  a  measurement,  larger  particles  may  pass 
further  out  and  still  trigger  a  measurement,  thus  defining  a  larger  active  area. 

5.2.1  On-Axis 

The  aperturing  effect  of  the  photodetector  pinhole  was  designed  in  this  system  to  occur 
outside  the  region  of  beam  crossover  so  that  visibilities  of  particles  being  sized  would  not  be 
altered  by  the  pinhole.  The  size  of  the  cross-sectional  area  in  this  case  is  determined , 
primarily  by  the  fall-off  in  illuminating  intensity  contained  in  the  exponential  factor  in  Eq./ 
(1)  and  can  be  calculated  as  follows. 


18 


AEDC-TR-8Q-50 


The  a-c  magnitude  signal,  separated  out  in  the  electronics,  is  used  as  the  reference  signal 
to  be  compared  with  the  selected  threshold  to  trigger  a  measurement.  The  a-c  (cosine) 
magnitude  is  here  designated  ia  and  from  Eq.  (1)  is  given  by 

=  Vo  ‘po  (D>  exP  [~2/bo  (*2  +  y2  +  z2  “2/4)]  (11) 

where  Voip^D)  is  the  a-c  magnitude  for  a  centered  particle  of  diameter  D  and  is  given  by  the 
curve  in  Fig.  4.  For  particles  moving  in  the  y-direction  (defined  in  Fig.  1),  ia  reaches  a 
maximum  value  where  the  particle  crosses  the  xz  plane  (y  =  0),  and  the  active  cross-sectional 
area  is  determined  in  this  plane. 

Following  a  derivation  in  Ref.  (12),  in  the  y  =  0  plane  Eq.  (11)  may  be  written 

i,  -  v0  ipo  (D)  exp  j-2  [(,>„’)-  (*2/b !)]}  (12) 

y 

where  b3  =  2b0/a.  This  equation  yields  elliptical  contours  of  constant  signal  amplitude,  ia, 

i 

in  the  y  =  0  plane,  described  by 

(-Ik)2  *  («A»)S  ■ (n  [(v°  ‘po <°>/‘.)«]  <13> 

and  the  area  within  the  contour  specified  by  setting  ia  equal  to  the  threshold  current,  io, 
gives  the  cross-sectional  area.  A,  for  a. particle  of  diameter  D: 

A  .  „b„l>al»  [(v.S.0»/.e)*]  (14) 


5.2.2  Off-Axis 

In  the  off-axis  system,  the  aperturing  effect  of  the  photodetector  slit  is  the  primary 
mechanism  which  determines  the  z-direction  dimension  of  the  active  cross-sectional  area,  A. 
The  fall-off  in  collected  light  with  the  coordinate  z  for  particles  passing  up  or  downstream  of 
the  z  =  0  point  is  shown  in  Fig.  5.  This  response  curve  is  approximately  Gaussian,  with 
standard  deviation  a  —  0.412  mm.  The  curve  has  fallen  to  1/e2  times  its  maximum  value  at  z 
—  2 a  =  0.824  mm,  much  faster  than  the  fall-off  in  illuminating  intensity  with  z  given  by  Eq. 
(1).  The  containment  of  the  active  cross-sectional  area  to  a  region  near  the  center  of  the 
beam  crossover  region  means  that  the  maximum  pedestal  ip  for  particles  moving  in  the  y 
direction  occurs,  as  does  the  maximum  a-c  component,  ia,  at  y  =  0.  The  Gaussian  shape  of 
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the  response  curve  of  Fig.  5  further  leads  to  the  approximation  that  the  pedestal  ip  in  the  y  = 
0  plane  may  be  given  by 


i=i 


pa 


(D)  expj-2  [(x2/b2)-(z2/b2) 


(15) 


which  is  similar  to  Eq.  (L2)  for  the  a-c  magnitude,  ia,  in  the  on-axis  system,  but  where  bj  is 
here  equal  to  the  experimental  value  0.824  mm. 


As  in  Eq.  (14)  above,  the  cross-sectional  area,  A,  is  given  by 

A  -  'b„M"[(<Po<“>/’ic)*]  06) 

where  the  pedestal  current  is  used  as  the  reference  signal  to  be  compared  with  the  threshold 
current,  iG. 

A  plot  of  the  cross-section,  A,  versus  diameter,  D,  calculated  using  Eq.  (16)  is  shown  as 
the  solid  line  in  Fig.  8.  The  curve  has  been  smoothed  below  10  pirn  to  eliminate  oscillations 
that  would  otherwise  occur  due  to  the  oscillations  in  the  ipo  versus  D  curve  as  shown  in  Fig, 
6.  The  system’s  resolution  below  10  /im  is  not  sufficient  to  make  accounting  for  such 
oscillations  necessary. 


Figure  8,  Sample  volume  cross-sectional  area  versus  particle 
diameter  for  off-axis  system. 
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If,  in  addition  to  the  lower  limit  on  the  pedestal  io.  an  upper  limit  is  also  placed,  say  ipmx> 
such  that  signals  whose  pedestals  exceed  ipra*  are  rejected,  then  the  cross-section.  A,  becomes 
the  annular  area  between  two  ellipses;  the  area  within  an  inner  ellipse  must  be  subtracted 
since  particles  passing  there  are  unacceptable  because  of  too  great  a  signal  magnitude.  The 
area.  A,  becomes 

A  =  n  bo  b3  £njjjpo  “  *  bo  b3  En  |[lpo  (D)/lPmx  } 

=  -b0  b32n  [(ipmx/iG)"*]  (I?) 

The  area  A  becomes  a  constant  given  by  Eq.  (17)  for  droplet  diameters  such  that 

ipo (D)  -  'pmx 

This  case  is  shown  as  the  dotted  line  in  Fig.  8.  Although  this  result  of  a  constant  sample 
volume  independent  of  particle  size  is  based  on  a  mathematical  model,  it  should  hold  for  a 
useful  range  of  particle  diameters.  It  should  also  be  noted  that  the  sensitive  cross  section,  A, 
is  annular,  which  means  that  large  particles  passing  through  the  center  of  the  crossover 
region  are  rejected  for  measurement,  although  this  region  is  normally  where  the  best 
measurements  relating  visibility  to  size  are  made. 

'  5.3  CORRECTIONS  SUMMARY 

In  summary,  the  correction  to  each  histogram  raw  bin  count  is  made  as  follows: 

( -  - 1 -  particles  /cm3Axm  (18) 

Vad^  A.v.^AD).  ‘ 


where 

( is  the  corrected  bin  count  (i.e.,  the  ith  bar  height), 

Vad;. 

Nj  is  the  raw  bin  count  for  the  ith  bin, 

A.  is  the  sample  volume  cross-sectional  area  in  square  centimeters  for  the 
particle  size  corresponding  to  the  ith  bin, 

Vj  is  the  velocity  in  centimeters  per  second  of  those  particles  whose  size 
corresponds  to  the  ith  bin, 
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tr  is  the  time  during  which  the  signal  processor  is  online,  sampling  the  flow, 
and 


( AD  )j  is  the  width  of  the  ith  bin  in  micrometers. 

If  only  the  shape  of  the  distribution  is  desired  (i.e.,  only  the  relative  frequency  of 
occurrence  of  each  size),  rather  than  an  absolute  count  in  particles/cm3,  and  if  it  is  assumed 
that  all  particles  have  the  same  velocity,  then  the  factor  Vjtr  in  Eq.  (18)  may  be  omitted. 
Further,  if  bin  widths  are  all  the  same  on  the  size  axis  and  the  cross-sectional  area.  A,  is  also 
a  constant,  then  no  corrections  are  required  to  obtain  a  relative  distribution  curve.  In  fact, 
for  the  off-axis  system  being  described  in  this  report,  these  factors  are  approximately 
constant,  and  the  raw  data  may  be  presented  alone  to  give  relative  distribution  curves. 

6.0  MONOSIZE  RESPONSE:  THEORETICAL  CONSIDERATIONS 

An  important  characteristic  of  a  particle-sizing  system  is  its  response  to  a  monosize 
(monodisperse)  flow  of  particles.  The  narrowness  of  this  response  determines  system 
resolution  (i.e.,  the  ability  of  the  system  to  distinguish  a  bimodal  distribution  when  the  two 
modes,  or  peaks,  are  closely  spaced).  Correspondingly,  it  determines  the  capability  of  the 
system  to  reproduce  a  steeply  rising  or  falling  distribution  curve. 

The  Berglund-Liu  vibrating  orifice  droplet  generator  is  capable  of  producing  a  highly 
monosize  flow  of  droplets  down  to  diameters  on  the  order  of  1  fim;  this  capability  enables 
one  to  conveniently  observe  the  monosize  response  while  calibrating  the  system.  When  the 
monosize  response  is  measured,  it  is  important  that  the  test  particles  be  permitted  to  flow 
freely  throughout  the  active  measurement  region  of  the  instrument,  rather  than  being 
constrained  to  flow  only  through  some  one  favorable  point. 

When  an  actual  distribution,  f(D),  is  measured,  the  distribution  recorded  by  the 
instrument  (i.e.,  the  system  output),  say  g(D),  is  given  by  the  convolution  of  the  actual 
distribution,  f(D),  and  the  monosize  response  of  the  system,  h(D): 

6(11)  =  /*  h(D-D0)f(Do)dDo  (19) 

In  this  formulation,  h(D)  is  the  response  to  a  “unit”  monosize  flow;  i.e.,  there  is  one 

no 

particle/cm3  in  the  flow:  J  h(D)dD  =  1.  It  is  also  assumed  that  the  monosize  response  is 

-  oc 

independent  of  the  diameter  of  the  monosize  droplets,  although  this  is  probably  not  actually 
the  case. 
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An  example  illustrating  Eq.  (19)  is  shown  in  Figs.  9a  and  b.  A  hypothetical  instrument 
monosize  response  is  shown  in  Fig.  9a,  here  taken  simply  to  be  a  Gaussian  response.  A 
hypothetical  log-normal  size  distribution,  f(D),  is  graphed  as  the  solid  curve  in  Fig.  9b.  The 
result  of  numerically  convolving  f(D)  and  h(D)  is  the  dashed  curve,  g(D),  shown  in  Fig.  9b. 
It  is  seen  that  the  instrument  in  this  case  can  reproduce  quite  well  the  right-hand  slope  of 
f(D)  but  cannot  reproduce  the  steeper,  left-hand  slope.  The  monosize  response  curves  for 
the  instruments  described  in  this  report  are  discussed  in  Section  7.1. 


D0,  pm  D0.  pm 

a.  Instrument  response  to  a  monosize  b.  Instrument  response,  g(D),  given  by 
droplet  flow  convolution  of  monosize  response 

and  size  distribution,  f(D) 

Figure  9.  Instrument  response. 

7.0  DATA  RESULTS 

7.1  LABORATORY  DATA 

Visibility  measurements  with  the  on-axis  system  were  obtained  from  glass  spheres 
(microbeads)  placed  on  a  glass  slide,  and  from  water  droplets  and  oleic  acid  droplets 
produced  by  a  Berglund-Liu  monosize  droplet  generator.  Experimental  data  shown  in  Fig.  2 
are  examples  of  glass  bead  data. 

The  data  acquisition  electronics  used  with  the  on-axis  system  measured  the  visibilities, 
converted  to  particle  size  through  the  solid  curve  of  Fig.  2,  and  then  displayed  a  size 


23 


AEDC-TR-80-S0 


distribution  histogram.  Figure  10  shows  the  result  of  measuring  100  droplets  each  of  five 
different  sizes  produced  by  the  droplet  generator.  The  drops  flowed  freely  throughout  the 
beam  crossover  region  during  the  test  and  were  not  constrained  to  a  single,  favorable  path. 
In  effect.  Fig.  10  shows  the  monosize  response  (see  Section  6,0)  of  the  instrument  at  five 
different  pdrticle  diameters.  The  monosize  response  is  quite  narrow  for  this  on-axis  system, 
giving  good  resolution,  adequate  to  easily  resolve  sizes  5  /*m  apart  in  this  range. 

Droplet 

Diameters,  pm 


Figure  10.  On-axis  system  response  to  five  different  monosiza 
droplet  flows. 


The  experimental  data  of  Fig.  3  are  examples  of  laboratory  data  taken  using  the  off-axis 
system  and  water  and  oleic  acid  droplets  produced  by  the  Bergiund-Liu  generator. 

The  monosize  response  of  the  off-axis  system  was  broader  than  that  of  the  on-axis 
system.  An  example  is  shown  in  Fig.  11,  for  a  monosize  stream  of  droplet  diameter  =  9.2 
/un.  The  broader  response  was  anticipated  since  the  off-axis  system  achieves  its  smaller 
observation  volume  by  the  photodetector  slit  aperture’s  partially  or  completely  blocking 


24 


AEDC-TB -80-50 


scattered  light  originating  anywhere  except  from  a  small  region  at  the  center  of  the  beam 
crossover  region.  The  resultant  partial  blocking  of  the  scattered  light  from  many  signals 
alters  the  visibility  for  these  signals  from  the  calibration  curve  of  Fig.  3  and  gives  a 
broadened  monosize  response. 


D,  pm 


Figure  11.  Off-axis  system  response  to  monosize  droplet  flow,  ■ 
diameter  =  9.2  jum. 


7.2  WIND  TUNNEL  DATA 

The  on-axis  system  was  installed  first  for  testing  in  the  ETF  RID  icing  studies  wind 
tunnel  to  obtain  water  droplet  measurements.  The  tunnel  has  a  3-ft  outside  diameter  with 
the  droplets  carried  within  a  12-in. -diam  flow.  Droplet  diameters  were  measured  in  the  range 
from  5  to  80  /im.  Particle  number  densities  ranged  up  to  a  few  hundred  particles/cm3. 
Velocities  in  the  flow  exceeded  100  m/sec.  Measurements  produced  inconsistent  results, 
which  are  attributed  to  a  high  probability  that  there  was  more  than  one  particle  at  a  time  in 
the  observation  volume. 
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The  off-axis  system  was  designed  to  provide  the  necessary  smaller  observation  volume 
and  was  installed  in  the  tunnel  after  laboratory  calibration.  An  example  of  a  data  histogram 
obtained  for  a  particular  set  of  tunnel  and  spray  nozzle  conditions  is  shown  in  Fig.  12.  The 
histogram  shown  displays  raw  data.  Data  corrections  discussed  in  Section  5.0  of  this  report 
were  not  necessary  to  show  the  relative  shape  of  the  distribution,  since  bin  widths  are 
uniform  along  the  size  axis  (due  to  the  linear  calibration  curve  shown  in  Fig.  3)  and  since  the 
sample  volume  cross-sectional  area  is  approximately  constant  for  sizes  greater  than  10  /xm 
(shown  as  the  dashed  line  in  Fig.  8).  To  attempt  to  correct  for  cross-sectional  area  variations 
below  10  f/m  would  only  have  served  to  disproportionately  amplify  noise  counts  in  this  case. 


Figure  12.  Log-normal  curve  fit  to  measured  tunnel 
size  distribution. 


The  sizes  of  spray  droplets  from  a  nozzle  are  often  modeled  by  a  log-normal  distribution, 
given  by 


/  (n  D  -  V.n  a\  \ 

\  f"  «2  /  . 


An 

AD 


exp 
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where  n  is  the  particle  number  density,  D  is  the  particle  diameter,  aj  is  the  geometric  mean, 
and  az  is  the  geometric  standard  deviation.  The  solid  curve  in  Fig.  12  is  a  log-normal  curve 
with  parameters  chosen  to  fit  the  experimental  data.  Parameters  used  were  lnai  =  2.4  and 
In  az  =  0.4.  A  multiplying  factor  was  then  used  to  give  the  log-normal  curve  the  same 
height  as  the  histogram. 

The  fit  fails  above  about  30  (im,  possibly  because  of  droplets  combining  to  form  bigger 
drops  after  an  initial  log-normal  distribution  was  produced  by  the  nozzle.  The  fit  also  fails 
for  droplets  smaller  than  about  6  jim,  which  would  be  expected  for  a  system  having  a 
monosize  response  as  shown  in  Fig.  11,  where  system  resolution  is  not  sufficient  to  follow 
the  steep  rise  of  the  log-normal  model.  (See  the  example  calculation  illustrated  in  Fig.  9.) 

An  alternative  model  sometimes  used  for  water  sprays  is  the  Nukiyama-Tanasawa 
distribution,  given  by 

AiL  =  a  D2  exp  bD^l  (21) 

AD  L  J 

t 

where  a,  b,  and  q  are  parameters  chosen  to  fit  the  experimental  data.  The  solid  curve  in  Fig. 
13  shows  a  Nukiyama-Tanasawa  curve  fit  to  the  same  data  histogram  as  that  shown  in  Fig. 
12.  Parameters  used  were  b  =  6.67  x  Id4  and  q  =  3,  with  the  factor  “a”  value  chosen  to 
give  the  curve  the  same  height  as  the  histogram.  The  histogram  seems  better  fit  by  the  log¬ 
normal  model  than  by  the  Nukiyama-Tanasawa  in  this  case. 

The  data  agree  reasonably  well  with  size  distributions  determined  from  holograms  taken 
of  the  droplet  flow  under  similar  tunnel  and  nozzle  conditions.  For  comparison  with  the 
holographic  data,  the  PS1  histogram  bin  widths  were  changed  from  1.5  to  4.5  fim;  i.e.,  raw 
counts  from  every  three  bins  were  combined  (summed)  into  one  equivalent  bin  having  a 
width  equal  to  the  combination  of  three  original  bins.  This  provides  better  comparison  with 
the  holographic  data  histograms,  which  have  bins  5  jim  wide. 

Figures  14  and  15  each  show  PSI  data  as  a  solid  line  and  holographic  data  as  a  broken 
line.  Each  holographic  histogram  has  been  multiplied  by  a  factor  to  approximately  equalize 
its  height  with  the  corresponding  PSI  histogram,  since  it  is  the  relative  distribution  that  is 
being  compared.  In  Fig.  14  the  PSI  data  were  taken  for  a  water  flow  of  0.1  gal/min  and  a 
nozzle  air  pressure  of  29  psi.  For  the  holographic  data,  the  conditions  were  similar,  being 
0.111  gal/min  and  50.31  psi. 

In  Fig.  15  the  conditions  for  the  PSI  data  were  0.044  gal/min  and  144.09  psi  and  for  the 
holographic  data  were  0.047  gal/min  and  112.5  psi. 
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Figure  13.  Nukiyama-Tanasawa  curve  fit  to  measured  tunnel 
size  distribution. 


Diameter,  pm 

Figure  14.  Comparison  of  particle  size 
interferometer  with 
holography  data. 
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The  comparison  of  PSI  and  holographic  data  is  reasonably  good.  It  is  assumed  that  the 
reason  the  PSI  shows  a  higher  count  of  particles  less  than  10  pm  in  diameter  is  that 
holography  is  less  sensitive  to  the  smaller  particles  in  that  range. 


Figure  15.  Comparison  of  particle  size  interferometer  with 
holography  data. 

8.0  CONCLUSIONS 

A  particle-sizing  interferometer  was  designed  and  built  for  use  in  the  AEDC  subscale 
icing  studies  facility,  RID.  Two  designs  were  evaluated,  one  with  forward-scatter,  on-axis 
light  collection,  and  one  with  collection  of  the  light  from  12  deg  off-axis  in  the  forward 
direction.  Laboratory  evaluation  showed  the  on-axis  system  to  have  the  better  capability  to 
resolve  between  closely  spaced  sizes.  The  requirement  of  a  small  observation  volume  to 
discriminate  against  multiple-particle  signals,  however,  dictated  the  use  of  the  off-axis 
system  to  obtain  correct  results  in  the  wind  tunnel.  Data  obtained  in  the  tunnel  compared 
favorably  with  holographic  data  obtained  under  similar  spray  nozzle  conditions  of  air 
pressure  and  water  flow. 

Water  droplets  were  sized  in  the  range  from  5  to  80  fim  while  droplet  velocities  measured 
in  the  air  flow  were  in  excess  of  100  m/sec.  Number  densities  in  the  flow  were  a  few  hundred 
droplets  per  cubic  centimeter. 

Particle  size  interferometry  offers  an  attractive  means  of  in  situ  sizing  of  droplets  in  the 
range  of  diameters  from  a  few  micrometers  to  a  few  hundred  micrometers.  The  method 
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features  rapid,  online  data  reduction  and  simultaneous  velocity  measurement. However,  it 
shares  with  other  single-particle,  in  situ  instruments  the  requirement  for  low  to  moderate 
particle  number  densities. 

Further  investigation  will  be  required  to  define  the  applicability  of  particle  sizing 
interferometry  to  irregularly  shaped  particles. 
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